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Abstract. We present a search for Trojan companions to 25 transiting exoplanets. We use the 
technique of Ford & Gaudi, in which a difference is sought between the observed transit time and 
the transit time that is calculated by fitting a two-body Keplerian orbit to the radial-velocity 
data. This technique is sensitive to the imbalance of mass at the L4/L5 points of the planet- 
star orbit. No companions were detected. The median 2cr upper limit is 60 M®, and the most 
constraining limit is 2.5 for the case of GJ 436. 
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1. Introduction 

Trojan companions are bodies in a 1:1 mean-motion resonance with a planet, librating around 
one of the two Lagrange points (L4 and L5) of the planet's orbit around the star. Several methods 
have been proposed to detect Trojan companions to exoplanets (Laughlin & Chambers 2002, 
CroU et al. 2007, Ford & Holman 2007, Ford & Gaudi 2006). In this paper, we present a search 
for Trojan companions to 25 known transiting exoplanets for which suitable data are available, 
using the method of Ford & Gaudi (2006, hereafter, "FG"). An important virtue of this method 
is that a sensitive search for Trojans can be performed using only the RV and photometric data 
that are routinely obtained while confirming transit candidates and characterizing the planets. 

2. Method 

The basic idea of the FG method is to compare the photometrically observed transit time 
{to) with the expected transit time {tc) that is calculated by fitting a two-body Keplerian orbit 
to the RV data, i.e. assuming no Trojan. The presence of a Trojan companion as a third body 
would cause a timing offset At = to — tc- 

In the case of a circular orbit, if there is no Trojan companion, the force vector on the star 
points directly at the planet and the observed transit time to coincides with the time tv when 
the orbital velocity of the star is in the plane of the sky (i.e., the time of RV null). If instead there 
is a single Trojan located at the L4 or L5 Lagrange point (or librating with a small amplitude), 
then the force vector on the star does not point directly at the planet; it is displaced in angle 
toward the Trojan companion. As a result, to occurs earlier or later than ty- The magnitude of 
to —tv is proportional to the Trojan mass ttit, for small values of the Trojan-to-planet mass ratio 
(Ford & Gaudi 2006): M ~ ±37.5 min (P/3 days)(mT/10 Me)(0.5 Mj/(mp + mr)) (1). The 
positive sign corresponds to a mass excess at the L4 point (leading the planet) while the negative 
sign corresponds to that at the L5 point (lagging the planet). For an eccentric two-body orbit, 
tc does not generally coincide with tv- We calculate tc by fitting a two-body Keplerian orbit to 
the RV data and calculating the expected transit time based on the fitted orbital parameters. 

3. Data Analysis 

The RV data were taken from the available literature on each system. We fitted the Keplerian 
model to the RV data using a Markov Chain Monte Carlo (MCMC) technique, employing a 
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Metropolis-Hastings algorithm (see, e.g., Ford 2005). In general, uniform priors were used for all 
parameters. A prior constraint on ecosoj was used where secondary eclipse measurements are 
available in literature. A single chain of ~ 10® links was used for each system. For each parameter, 
we found the mode of the a posteriori distribution, and the 68.3% confidence interval, defined 
as the range that excludes 15.9% of the probability at each extreme of the distribution. We 
thus obtained tc from fitting the RV data. We determined to using the most precise published 
photometric ephemeris, and then computed At = to — tc- The results for At are translated 
into constraints on the Trojan mass rriT using Eq. (1) for circular orbits, and using equivalent 
expressions for eccentric orbits obtained by numerical integrations of three-body orbits. 



4. Results 

The 95.4% (2cr) upper limits on tut and on mT/mp are shown in Figure 1. The median upper 
limit on mx is 60 M®, with the most constraining limit of 2.5 M® holding for GJ 436. This 
powerful upper limit is possible for this case because of the small stellar and planetary masses, 
and the copious RV data available for this system. The median upper limit on the mass ratio 
niT/mp is 0.1, with a best-case value of 0.015 for HD 17156. The powerful constraint in this 
case is a result of the unusually long orbital period of 21 days and plentiful precise RV data. 
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Figure 1. 95.4 % upper-limits on trojans masses and on trojan-planet mass ratios. The systems 
are ordered from least constrained to most well constrained, going from left to right of the figure. 

In all cases but two, the result for At was consistent with zero within 2a. The exceptions 
were CoRoT-Exo-2 and WASP-2, for which M = SOl" and -142t|!^ minutes respectively. The 
results for these two systems are worth following up with additional RV data. However, in a 
sample of 25 systems, even if At is always consistent with zero, one expects approximately one 
2(7 outlier. There were no cases in which At was inconsistent with zero at the 3cr level. Hence, 
we conclude that there is no compelling evidence for a Trojan companion in this ensemble. 
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